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Abstract 

We consider a plane electromagnetic wave incident on a periodic stack of dielectric layers. One 
of the alternating layers has an anisotropic refractive index with an oblique orientation of the 
principal axis relative to the normal to the layers. It was shown recently that an obliquely incident 
light, upon entering such a periodic stack, can be converted into an abnormal axially frozen mode 
with drastically enhanced amplitude and zero normal component of the group velocity. The stack 
reflectivity at this point can be very low, implying nearly total conversion of the incident light into 
the frozen mode with huge energy density, compared to that of the incident light. Supposedly, 
the frozen mode regime requires strong birefringence in the anisotropic layers - by an order of 
magnitude stronger than that available in common anisotropic dielectric materials. In this paper 
we show how to overcome the above problem by exploiting higher frequency bands of the photonic 
spectrum. We prove that a robust frozen mode regime at optical wavelengths can be realized in 
stacks composed of common anisotropic materials, such as YVO4, LiNb, CaCO^, and the like. 
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I. INTRODUCTION 



In photonic crystals, the speed of light is defined as the wave group velocity, u, 



u = du/dk, 



(1) 



where k is the Bloch wave vector and uj = to \ kj is the respective frequency. At certain 
frequencies, the dispersion relation uj ( k ) of a photonic crystal develops stationary points 



in the vicinity of which the group velocity vanishes. Zero group velocity usually implies 
that the respective Bloch eigenmode does not transfer electromagnetic energy. Indeed, with 



where W is the electromagnetic energy density, associated with this mode. If W is limited, 
then the group velocity u and the energy flux S vanish simultaneously at any stationary 
point (J2J) of the dispersion relation. Such modes are commonly referred to as slow modes, 
or slow light. Examples of different stationary points (J2J) are shown in Fig. 0, where each 
of the respective frequencies u a , u>b, uj 9 and c^o is associated with slow light. 

A common problem with slow modes is that most of them cannot be excited in semi- 
infinite photonic crystal by incident light. Indeed, consider plane monochromatic wave 
incident on a semi-infinite photonic crystal with the electromagnetic dispersion relation 
shown in Fig. ^ If the frequency uj is close to the band edge frequency uj 9 in Fig. ^ then 
the incident wave will be totally reflected back into space, as seen in Fig. |21 

In another case, where the incident wave frequency is close to the characteristic value 
uj a or Ub in Fig. Q some portion of the incident wave will be transmitted in the photonic 
crystal, but none in the form of the slow mode corresponding to the respective stationary 
point. This means, for example, that at the frequency uj a , the transmitted light is a Bloch 
wave with the finite group velocity and the wave number different from that corresponding 
to the point a in Fig. ^ 

Let us turn now to the stationary inflection point of the dispersion relation in Fig. d 
where both the first and the second derivatives of the frequency uj with respect to k vanish, 




duj/dk = 



(2) 



certain reservations, the energy flux S of a Bloch mode is 



S = Wu 



(3) 



2 




wave number k 

FIG. 1: An example of electromagnetic dispersion relation u(k) with various stationary points: 
(i) extreme points a and b of the respective spectral branches, (ii) a photonic band edge g, (hi) a 
stationary inflection point 0. Each stationary point is associated with slow light. The wave number 
k and the frequency u are given in units of 1/L and c/L, respectively. L is the primitive translation 
of the periodic array. 



while the third derivative is finite 



. , , du 
at u> = uj and k = ko : — — 

ok 



0; 



d 2 u 



(4) 



dk 2 dk 3 

The frequency ujq of stationary inflection point is associated with the so-called frozen mode 
regime |2|, y] . In such a case, the incident plane wave can be transmitted into the photonic 
crystal with little reflection, as seen in Fig. El Having entered the photonic slab, the 
light is 100% converted into the slow mode with infinitesimal group velocity and drastically 
enhanced amplitude. Under the frozen mode regime, vanishingly small group velocity u in 
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FIG. 2: Transmittance r as a function of incident light frequency u for the semi-infinite photonic 
slab with the dispersion relation presented in Fig. 1. The characteristic frequencies uj a , ^b, ^o, and 
ujg are associated with the respective stationary points in Fig. 1. At u > uj g (within the photonic 
band gap) the incident light is totally reflected by the slab. 

Eq. (J2J) is offset by a huge value of the energy density W. As the result, the energy flux 
(j2J) associated with the frozen mode remains finite and comparable with that of the incident 
wave. In the vicinity of the frozen mode frequency u , the electromagnetic energy density 
W associated with the slow (frozen) mode displays a resonance-like behavior 

^-HK')" /3 (---or 2/3 , (5) 

where Sj is the fixed energy flux of the incident wave, r is the portion of the incident light 
transmitted into the semi-infinite photonic crystal, and 
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Remarkably, the transmittance r at 00 ~ c^o remains finite and may even be close to unity, as 
shown in Fig. El The latter implies that a significant portion of the incident light is converted 
into the frozen mode with nearly zero group velocity and huge amplitude, compared to that 
of the incident wave. In reality, the electromagnetic energy density W of the frozen mode 
will be limited by such factors as absorption, nonlinear effects, imperfection of the periodic 
dielectric array, deviation of the incident radiation from a perfect plane monochromatic 
wave, finiteness of the photonic slab dimensions, etc. Still, with all these limitations in 
place, the frozen mode regime can be very attractive for various practical applications. 

From now on we restrict ourselves to the case of lossless periodic layered media (periodic 
stacks), which can be viewed as one dimensional photonic crystals. According to |2j, at 
normal incidence, the frozen mode regime in a periodic stack can only occur if some of the 
layers display sufficiently strong circular birefringence (Faraday rotation). In addition, each 
unit cell of the periodic layered array must contain at least two layers with significant and 
misaligned in-plane anisotropy. If the above conditions are not met, the electromagnetic 
dispersion relation of the periodic stack cannot develop stationary inflection point (J3J) and, 
therefore, cannot support the frozen mode regime at normal incidence. At the microwave 
frequency range, one can find a number of materials meeting the above requirements. But 
at infrared and optical frequencies, the circular birefringence of known transparent magnetic 
materials becomes too small to support a robust frozen mode regime Q]. Since our prime 
interest here is with optics, we will explore the "non-magnetic" approach proposed in j^j]. 

According to |2j, the frozen mode regime can occur in nonmagnetic periodic stacks with 
special configurations requiring some layers to display appreciable oblique (neither in-plane, 
nor axial) anisotropy. On the other hand, at optical frequencies, all commercially available 
anisotropic dielectrics display substantially weaker anisotropy, compared to what would be 
the optimal value. According to [3j, too weak anisotropy can push the frozen mode frequency 
too close to the nearest band edge, resulting in almost total reflectance of the incident light. 
The high reflectance of the slab, in turn, implies very low efficiency of conversion of the 
incident light into the slow mode. In this paper we show that in fact, the negative effect of 
the weak anisotropy on the frozen mode regime can be completely overcome by proper design 
of the layered structure. As the result, a robust frozen mode regime at optical frequencies 
can be achieved in periodic stacks incorporating real anisotropic materials such as yttrium 
vanadate, lithium niobate, and the like, where the dielectric anisotropy is one-two orders of 
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magnitude short of the "optimal" value. The idea is to choose the parameters of the periodic 
stack so that a stationary inflection point associated with the frozen mode regime develops 
at higher photonic bands. For a given frequency range, this requires thicker dielectric layers, 
which could be an additional practical advantage. A side effect of using higher photonic 
bands is that the effective bandwidth of the frozen mode regime appears to be narrower, 
compared to the case of hypothetical materials with much stronger anisotropy used in j^] 
for numerical simulations. 

The practical development of frozen mode devices from such commodity materials could 
lead to revolutionary advances in optical computing, sensing, and information processing. 
When practically realized, such frozen mode structures would enable significant advances in 
all-optical information storage and processing (such as optical memory and buffer elements, 
optical delay lines) as well as optical sensing, lasing, and nonlinear optics. 

The paper is organized as follows. In the next section, we discuss in general terms the 
phenomenon of axially frozen mode. The detailed analysis of the mathematical aspects of the 
phenomenon can be found in j^J . Then, in section III, using a specific example of a periodic 
array incorporating yttrium vanadate, we demonstrate how a robust frozen mode regime 
at optical frequencies can be achieved in a practical setting involving weakly anisotropic 
materials. Finally, in Appendix, we briefly overview the electrodynamics of lossless layered 
media, introducing basic notations, definitions, and assumptions used in our computations. 



II. AXIALLY FROZEN MODE REGIME AT OBLIQUE INCIDENCE 

According to , in nonmagnetic periodic stacks, the simplest version of the frozen 
mode regime described in the introduction is impossible. Still, a more general phenomenon 
referred to as the axially frozen mode regime can occur. This section starts with a brief 
general discussion of the phenomenon. Then we turn to the particular case of a periodic 
stack incorporating yttrium vanadate layers. The reason we have chosen this particular 
material is because its optical properties is very similar to those of other common anisotropic 
dielectrics transparent at optical wavelength. 
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A. Basic definitions 



Consider a monochromatic plane wave obliquely incident on a periodic semi-infinite stack, 
as shown in Fig. El 

Semi-infinite 
periodic stack 
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FIG. 3: Light incident on a semi-infinite periodic layered medium. The arrows show the energy 
fluxes of the incident, reflected and transmitted waves, respectively. The transmitted wave is 
a superposition of two Bloch eigenmodes, each of which can be either propagating or evanescent. 
Only propagating modes can transfer electromagnetic energy in the z direction. 



Let \T/# and denote the incident, reflected and transmitted waves, respectively. 
Due to the boundary conditions pUj) . all three waves \T//, \T/# and \T/t must be assigned the 
same pair of tangential components k x , k y of the respective wave vector 0, ^| 



R 



1 T 



ki 



(6) 



y v / y v / y 

while their axial (normal) components k z can be different. Hereinafter, the symbol k z will 
refer only to the transmitted Bloch waves propagating inside the semi-infinite slab 



Inside periodic stack (at z > 0): k = {k x , k y , k z ) . 



(7) 



The value k z is defined up to a multiple of 2tt/L (the Brillouin zone), where L is the period 
of the layered structure. For given k x , k y and u, the value k z is found by solving the time- 
harmonic Maxwell equations (|37j) in periodic medium, that will be done in the following 
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sections. The result can be represented as the axial dispersion relation, which gives the 
correspondence between uj and k z at fixed k x , k y 

uj = uj (k z ) , at fixed k x , k y . (8) 

It can be more convenient to define the axial dispersion relation as the correspondence 
between uj and k z at fixed direction n of incident light propagation 

uj = uj (k z ) , at fixed n x , n y , (9) 



where the unit vector n can be expressed in terms of the tangential components © of the 
wave vector 




L 



FIG. 4: Periodic array of anisotropic dielectric layers (A) separated by gaps (B). The anisotropy 
axis of the dielectric material (the tetragonal axis, in the case of yttrium vanadate) makes an 
oblique angle with the z direction, normal to the layers. The stack parameters used in our numerical 
simulations are specified in formulae (30) and (31). 

The transmitted electromagnetic field inside the periodic layered medium is not a 
single Bloch mode, but it is a superposition of two Bloch modes with different polarizations 
and different values of k z . Of course, the tangential components k x , k y are the same for 
either transmitted Bloch mode and the incident wave, as stated by Eq. (JUJ). Generally, there 
are three possibilities (see the details in JjJ): 
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FIG. 5: The axial dispersion relation to(k z ) of the periodic stack in Fig. 4 at fixed direction n of 
light incidence (n x = n y = —0.493489). Small spectral asymmetry and small branch separation 
are due to weak anisotropy of yttrium vanadate. Both the asymmetry and the branch separation 
are more prononced in upper spectral bands. The wave number k z and the frequency u are given 
in units of 1/L and c/L, respectively. 

1. Both Bloch components of the transmitted wave are propagating modes, which 
means that the respective values of k z are real. For example, at u < u a and uo > ujf, in 
Fig. Elwe have two Bloch modes propagating inside the slab with two different group 
velocities u z > (double refraction). Note that propagating modes with u z < 0, as 
well as evanescent modes with S (k z ) < 0, do not contribute to the transmitted wave 
inside the semi-infinite stack in Fig. |H1 
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FIG. 6: A fragment of the axial dispersion relation uj (k z ) in Fig. 5, which includes the axial 
stationary inflection point at ujq = 9.164450223. 

2. Both Bloch components of are evanescent, which implies that the respective values 
of k z are complex with S (k z ) > 0. In particular, this is the case if the frequency uj 
falls into a photonic band gap (for example, at uj > uj 9 in Fig. [TJ. In such a case, the 
incident wave is totally reflected back to space. 

3. Of particular interest is the case when one of the Bloch components of the transmitted 
wave is a propagating mode (with u z > 0), while the other is an evanescent mode 
(with 3 (k z ) > 0) 

m T = ty pr + * e „. (11) 
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This is the case at the frequency range 

uj a < uj < Ub (12) 

in Fig. |U1 As the distance z from the slab/vacuum interface increases, the evanescent 
contribution ty ev decays as exp (— z*3 (k z )), and the resulting transmitted wave VP^ 
turns into a single propagating Bloch eigenmode ^ pr . 

Similarly to the case (J1J) of regular frozen mode, the axially frozen mode is associated 
with the axial stationary inflection point defined as 

at U = UJ ° : dk- z = 0; M = 0; M ^ a (13) 

The regular stationary inflection point is a particular case of (jlHj) . Example of axial 
dispersion relation displaying such a singularity is shown in Fig. In the vicinity of 
uq in Eq. (JT3|l . the electromagnetic field inside the slab is a superposition (fTTj) of 
one propagating and one evanescent Bloch components. As the frequency uj approaches 
the critical point (|T3j) . both contributions grow sharply, while remaining nearly equal and 
opposite in sign near the slab boundary (at z = 0), as illustrated in Fig. [7| 



a) Propagating mode b) Evanescent mode c) Resulting field 




9.16 9.17 9.16 9.17 9.16 9.17 

Frequency co Frequency co Frequency co 

FIG. 7: Destructive interference of the propagating and evanescent contributions to the resulting 
field at the slab/vacuum interface under the frozen mode regime: a) amplitude |^ pr (0)| 2 of 
the propagating component, b) amplitude \^f ev (0)| 2 of the evanescent component, c) resulting field 
amplitude \^>t (0)| 2 . The incident light has unity amplitude and TE polarization. 

Due to the destructive interference at the slab boundary, the resulting electromagnetic 
field at z = is small enough to satisfy the boundary condition (jlrjj) . As the distance z from 
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the slab boundary increases, the evanescent component ty ev decays exponentially, while the 
amplitude of the propagating component ^ pr remains constant and large, as shown in Figs. 
H> and Els. 



a) n = n Q , co = co Q b) n = n , co * co c) n * n Q , co = co Q 
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FIG. 8: The Amplitude \ty (z)\ 2 of the resulting field as function of the distance z from the slab 
boundary: a) the exact point of the frozen mode regime, b) the frequency is deviated from loq by 
10~ 3 c/L, c) the direction of light incidence is deviated from no in l)33jl by 10" 4 . The distance z is 
given in units of L. The incident light has unity amplitude and TE polarization. 



B. Energy density and energy flux of the axially frozen mode 

Let Si, Sr and St be the energy fluxes of the incident, reflected and transmitted waves, 
respectively. Within the frequency range (|12jh which includes the critical point ()13j) . the 
transmitted wave ^>t is a superposition (jlljl of propagating and evanescent components. 
Only the propagating component ty pr is responsible for the axial energy flux \Stj ■ 

The axial energy flux can also be expressed in terms of the axial component u z of the 
propagating mode group velocity and the energy density Wq associated with ^> pr 

(s^j = W u z oc \^ P r\ 2 u z , (14) 

The quantity Wo in Eq. (JHJ) can be interpreted as the electromagnetic energy density 
far from the slab interface, where the electromagnetic field \&t reduces to its propagating 
component ^ pr . In the vicinity of the axial stationary inflection point 1)130 . the energy 
density W , associated with the axially frozen mode, diverges, while u z — - > 0. As a result, 
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the vanishingly small u z in Eq. (|14|) is offset by a very large value of Wq. The theoretical 
analysis of the next section shows that the axial energy flux \J>tj i n Eq. f)14j) . along with 
the slab transmittance r in ([54)1 . remain finite even at = co> , where the axial component 
of the group velocity vanishes 

'S T ) >0, ifw 2 = 0. (15) 



The energy conservation consideration allows to find the asymptotic frequency depen- 
dence of the amplitude |^/ pr | 2 of the axially frozen mode in the vicinity of the critical point 
([13)1 . Indeed, in the vicinity of ujq, the axial dispersion relation can be approximated by a 
cubic parabola 

1 .»> n. i. \3 mi ( d 3uJ 



u - uqk, -oj'q (k z - ho) , ^"=( — ) . (16) 



fc=fco 



The 2; component of the group velocity is 

duj 1 , ; , (r- 



V " (a. - k ) 2 » V K') 1/3 (« - ^o) 2/3 • (17) 



<% z 2 

The Eq. f|17|) together with (JT4"jl yield the following asymptotic expression for the energy 
density W associated with the frozen mode 



or, equivalently, 



Wo * 6^ ^" 1/S (W - ^° r2/3 ' (18) 



^0 - fe) K")" 1/3 (" - ^o)- 2/3 • (19) 



6 2/3 

where Sj is the fixed energy flux of the incident wave and r is the transmittance coefficient 
defined in ()54j) . Remarkably, the transmittance r of the semi- infinite slab remains finite in 
the vicinity of the frozen mode frequency ujq, as seen in Figs. E]andE3 This implies that the 
electromagnetic energy density Wq associated with the frozen mode, as well as its amplitude 



l^/pj.) 2 , diverge as u — > u . In Figs. H] and El such a behavior is illustrated for two different 
incident light polarizations. 

In reality, under the axial frozen mode regime, the field amplitude inside the slab is 
limited by various physical factors mentioned earlier in this paper. Still, with all these 
limitations in place, the normal energy flux (jSrJ remains finite and comparable with that 
of the incident wave. The latter implies that a substantial portion of the incident wave is 
converted into the axially frozen mode with drastically enhanced amplitude and nearly zero 
axial component u z of the group velocity. In many respects, the phenomenon of axial frozen 
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FIG. 9: At the top: the red solid line shows the amplitude |^ P r| 2 of the propagating Bloch 
component of the transmitted wave *$>t from Eq. (11). The dotted green line shows the amplitude 
\^f e v\ of the evanescent component. Both contributions to *$>t diverge at the frozen mode frequency 
ojQ. At z S> L, the evanescent component decays and |^ P r| 2 (the solid line) represents the resulting 
amplitude |Wy| of the transmitted wave. By contrast, at uj < u a and u> > oJb, the transmitted wave 
\Py is a superposition of two propagating components (double refraction). At the bottom: the slab 
transmittance r vs. frequency u. At u = loo, the transmittance remains finite. The characteristic 
frequencies u> a , ujq and LOb are explained in Fig. 6. The incident wave has TE polarization and 
unity amplitude. 
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FIG. 10: The same as in Fig. 9, but the incident wave polarization is now TH. 

mode is similar to its particular case, the regular frozen mode, associated with the regular 
stationary inflection point (@J). 



C. Physical conditions for the axially frozen mode regime in layered media 

The physical conditions under which a non-magnetic layered structure can support the 
(axial) frozen mode regime can be grouped in two categories. The first one comprises 
several symmetry restrictions. The second category includes some basic qualitative recom- 
mendations which would ensure the robustness of the frozen mode regime, provided that the 
symmetry conditions for the regime are met. In what follows we briefly reiterate those con- 



15 



ditions and then show how they apply to periodic stacks incorporating some real dielectric 
materials. 

1. Symmetry conditions 

There are two fundamental necessary conditions for the frozen mode regime. The first 
one is that the Bloch dispersion relation oj (jtj in the periodic layered medium must display 
the so-called axial spectral asymmetry 

u (k x , k y , k x ) ^ uj (k x , k y , -k s ) . (20) 

As shown in jsj], this condition is necessary for the existence of the axial stationary inflection 
point (J13)) in the electromagnetic dispersion relation of an arbitrary periodic layered medium. 

The second necessary condition is that for the given direction k of wave propagation, the 
Bloch eigenmodes with different polarizations must have the same symmetry. In the case 
of oblique propagation in periodic layered media, the latter condition implies that for the 
given k, the Bloch eigenmodes are neither TE nor TM 

% is neither TE nor TM. (21) 

The condition (J2UJ) imposes certain restrictions on (i) the point symmetry group G of the 
periodic layered array and (ii) on the direction k of the transmitted wave propagation inside 
the layered medium. While the condition (J2*T|) may impose some additional restriction on 
the direction of k. 

The restriction on the symmetry of the periodic stack following from the requirement 
(|20|) of the axial spectral asymmetry is 

m z <£ G and 2 Z <£ G. (22) 

where m z is the mirror plane parallel to the layers, 2 Z is the 2-fold rotation about the z axis. 
An immediate consequence of the criterion (|22jl is that least one of the alternating layers of 
the periodic stack must be an anisotropic dielectric material with nonzero e xz and/or e yz , 
where the z direction is normal to the layers 

e xz ^ and/or e yz ^ (23) 



16 



Otherwise, the operation 2 Z will be present in the symmetry group G of the periodic stack. 
The condition (j20j) also imposes a restriction on the direction of wave propagation. Specif- 



neither parallel, nor perpendicular to the z direction 

k x ^ and/or k y ^ 0. 

The latter condition implies that the frozen mode regime cannot occur at normal incidence, 
regardless of the periodic stack geometry and composition. While the condition ()23|) implies 
that at least one of the alternating layers must be cut at an oblique angle relative to the 
principle axes of its permittivity tensor. If either of the above two conditions is not satisfied, 



which rules out the possibility of stationary inflection point and the frozen mode regime. 

If all the above necessary conditions are met, then the (axial) frozen mode regime is, at 
least, not forbidden by symmetry. More details on the symmetry aspects of the frozen mode 

n 

regime can be found in [3], Section II and |2j, Sections I and II. 
2. Additional physical requirements 

In practice, as soon as the symmetry conditions are met, one can almost certainly achieve 
the (axial) frozen mode regime at any desirable frequency to within certain frequency range. 
The frequency range is determined by the layers thicknesses and the dielectric materials 
used, while a specific value of uj within the range can be selected by the direction n of the 
light incidence. The problem is that unless the physical parameters of the stack layers lie 
within a certain range, the effects associated with the frozen mode regime can be insignificant 
or even practically undetectable. The basic guiding principle in choosing appropriate layer 
materials is "moderation". As soon as the "moderation" principle is observed, one can 
almost certainly achieve the frozen mode regime at prescribed frequency by choosing the 
right direction of light incidence. Specifically, those "moderation" conditions include: 

1. It is desirable that the ratio 



ically, the Bloch wave vector k must be oblique to the stack layers, which means that k is 





(24) 



An 



(25) 



n 
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of the birefringence and the refractive index of the material of the anisotropic layers 
lies somewhere between 2 and 10. If the anisotropy is extremely strong or too weak, 
the Bloch waves with different polarizations become virtually separated, which ex- 
cludes the possibility of a robust frozen mode regime. For example, in the practically 
important case of extremely small anisotropy, the two transmitted Bloch waves can 
be approximately classified as TE and TM modes, which is incompatible with the 
symmetry condition (|2*Tjl for the frozen mode regime. 

2. The dielectric contrast of the adjacent layers A and B should be significant, but not 
extreme. The ratio ua/ub anywhere between 1.5 and 20 would be appropriate. In 
addition, the dielectric contrast between the layers should match the ratio (|25j) in the 
anisotropic layers: weaker anisotropy would require weaker dielectric contrast between 
the layers. 

3. Typical layer thickness should be of the order of A/4n, where A = 2k/oj is the light 
wavelength in vacuum and n is the respective refractive index. In reality, the acceptable 
layer thickness can differ from X/4n by several times either way. But too thick layers 
would push the stationary inflection point to high-order frequency bands, while too 
thin layers would exclude the possibility of the frozen mode regime at the prescribed 
frequency range. 

The biggest challenge at optical frequencies poses the first condition, because most of the 
commercially available optical anisotropic crystals have the ratio An/n of only about 0.1. 
According to [sj], this would push the axial stationary inflection point (fT3*|) very close to the 
photonic band edge and make the photonic crystal almost 100% reflective. This indeed would 
be the case if we tried to realize the frozen mode regime at the lowest frequency band. But, in 
the next section we show that one can successfully solve this problem by moving to a higher 
frequency band. So, a robust axially frozen mode regime with almost complete conversion 
of the incident light into the frozen mode can be achieved with the commercially available 
anisotropic dielectric materials displaying the birefringence ratio (|23|) by one-two orders of 
magnitude smaller, compared to the optimal value used in for numerical simulations. 
The drawback though is that using higher photonic frequency bands narrows the bandwidth 
of the frozen mode regime by roughly an order of magnitude. 
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III. PERIODIC STACK INCORPORATING YTTRIUM VANADATE LAYERS 



The simplest non-magnetic periodic stack satisfying the symmetry conditions (|22|) and 
(J23j) and, therefore, capable of supporting the axial frozen mode regime, is shown in Fig. 
El It is composed of anisotropic A layers separated by empty gaps. The anisotropic dielec- 
tric layers must display nonzero components e xz and/or e yz . The simplest choice for the 
respective permittivity tensor is 



&xx £ X z 

Eyy 

E-xz £ zz 



(26) 



where the z axis coincides with the normal to the layers. 



Yt 



nm js] 



;rium vanadate is a tetragonal dielectric with the permittivity tensor, £yy, at A = 1550 



Eyy 



en 
e 22 
£ 33 



4.62 
4.62 
3.78 



(27) 



where the Cartesian axis X3 is chosen parallel to the crystallographic axis C4. In order to 
achieve a non-zero e xz component one has to rotate the tensor (|2*7jl about the y axis by an 
angle 9, different from and tt/2 9]. The result of the rotation is 



en cos 2 6 + £33 sin 2 9 (eu — £ 33 ) cos 9 sin 9 

E 22 

(exi — ^33) cos 9 sin 9 e 33 cos 2 9 + En sin 2 9 



(28) 



For numerical simulations we can choose, for instance, 



9 = tt/4. 

In this case, Eq. (l2*5|l together with (}2*Tj) yield 



(29) 



4.20 4.42 

4.62 
4.42 4.20 



(30) 



which is compatible with the required form (j26j). 
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Let d,A and ds denote the thickness of the A layers and the thickness of the gaps between 
them, respectively. For our numerical simulations we can choose 

d A = d B = L/2, (31) 

where L is the period of the layered array in Fig. HJ 

Let us reiterate that the parameters (|29|) and (|31|) are chosen at random. In practice, we 
can always adjust them so that the stack suits specific practical requirements. The struc- 
tural period L should be chosen so that the frozen mode regime occurs within a prescribed 
frequency range. Then, the direction of light incidence can be adjusted so that the axial 
frozen mode regime occurs exactly at a prescribed frequency ujq. 

Symmetry arguments similar to those presented in 3J] show that in the case of the periodic 
array in Fig. 0J the necessary conditions ()22j) and (}2*Tj) for the frozen mode regime are satisfied 
only if the direction n of the light incidence lies neither in the x — z, nor y — z plane 



n x 7^ and n y ^ 0. (32) 



A. Electromagnetic properties in the vicinity of axial frozen mode regime 

In Fig. El we presented the typical axial dispersion relation uj(k z ) of the periodic stack in 
Fig. H]at fixed direction n of light incidence. In Fig. El we chose 

n x = n y = n o, where n = —0.493489, (33) 

because this particular direction of light incidence produces the axially frozen mode regime 
at certain frequency ujq, shown in Fig. El Due to the relatively small anisotropy of yttrium 
vanadate, the axial dispersion relation in Fig. El displays rather weak asymmetry uj{k z ) ^ 
uj (—h z ), which makes it virtually impossible to develop a stationary inflection point at the 
lowest spectral branches. As we go to upper photonic bands, the situation improves. In Fig. 
El we present the enlarged fragment of the axial dispersion relation in Fig. El This fragment 
covers the boundary region between the forth and the fifth bands. At frequency 

c^ = 9.164450223c/L (34) 

one of the spectral branches develops axial stationary inflection point (|13|). associated with 
the possibility of the axial frozen mode regime. 
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Our numerical analysis based on the transfer matrix approach (the computational details 
are presented in the next section) indeed shows a very robust (axial) frozen mode regime in 
this setting, in spite of the fact that the dielectric anisotropy of yttrium vanadate is more 
than an order of magnitude short of the optimal value. The drawback though is that the 
frequency bandwidth of the effect is roughly an order of magnitude narrower, compared to 
what could be achieved with hypothetical materials displaying much stronger anisotropy at 
optical frequencies. 

Let us start with the results presented in Figs. M and [TUJ The bottom plots in both 
figures display the frequency dependence of the stack transmittance r for two different 
polarizations of incident light . Clearly, in the vicinity of the frozen mode frequency ujq, the 
transmittance remains significant, which implies that a significant portion of the incident 
radiation is converted into the frozen mode. The top plots in Figs. 01 and El display the 
amplitude of the two Bloch components of the transmitted wave ^t- The solid and dotted 
lines correspond to the propagating and evanescent Bloch components, respectively. In the 
vicinity of the frozen mode frequency (at uj a < uj < ujf,), there is one propagating component 
(\l/ pr ) and one evanescent component (^ ev ), each of which diverges as uj — > ujq, in accordance 
with Eq. (|19j). At uj < uj a and uj > uj^, the transmitted wave ^>t is a superposition of two 
propagating components with different group velocities, which constitutes the phenomenon 
of double refraction. The characteristic frequencies in Figs. El and ^] are explained in Fig. 



Fig. [3 shows the frequency dependence of the resulting field amplitude |\Pt(0)| at the 
slab boundary, along with the amplitudes of its propagating and evanescent components 



uj — > ujq, their superposition (0) remains finite, to meet the boundary conditions at 
z = 0. As we move further away from the slab boundary, the evanescent component dies out, 
while the propagating mode ^ w (z) remains constant and large. As a result, the destructive 
interference of ty pr (z) and ^> ev (z) is removed, and the resulting field amplitude |\1/t(-2)| 2 
grows and approaches the value |\I/pr| 2 - This scenario is illustrated in Fig. and Eh If the 
frequency uj of the incident wave and its direction of propagation n exactly correspond to 
the critical values ujq and n , then the electromagnetic field ty? i z ) inside the semi- infinite 
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Although each of the two Bloch contributions to (0) diverges as 



21 



stack is described by a linearly diverging non-Bloch Floquet eigenmode 

|^t| 2 oc z 2 , 

as shown in Fig. |H]a. 

By way of example, let us present the actual geometrical parameters of the stack sup- 
porting the axially frozen mode regime for the case of infrared light with A = 1550 nm and 
the direction of incidence (|33|) . The expression (J34j) for the frozen mode frequency yield 

d = L/2 = 1130.4 nm. 

In practice, we do not have to adjust the layer thicknesses in order to achieve the frozen 
mode regime at a prescribed wavelength. Instead, we can tune the system into the axially 
frozen mode regime by adjusting the direction n of light incidence. 
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IV. APPENDIX. SCATTERING PROBLEM FOR ANISOTROPIC SEMI- 
INFINITE STACK 

In this section we briefly discuss the standard procedure we use to do the electrodynamics 
of stratified media incorporating anisotropic layers. For more details, see, for example, 
P. liol 11, 3| and references therein. 



A. Time-harmonic Maxwell equations in periodic layered media 

Our consideration is based on time-harmonic Maxwell equations 

VxE(i, y, z) = i^B (x, y,z), V x H (x, y, z) = -i^D 0, V, z) (35) 

with linear constitutive relations 

D (x, y,z) =e (z) E (ar, y,z), B (x, y,z)=ft (z) H (x, y, z) . (36) 

In layered media, the tensors e and fi in Eq. (j3T?j) depend on a single Cartesian coordinate 
z. Plugging Eq. (JBTIj) into (|33j) yields 

V x E (x, y, z) = i-pL (z) H (x, y, z) , VxH (x, y, z) = -i-e (z) E (x, y, z) . (37) 

c c 

Solutions for Eq. (j3~Tj) are sought in the following form 

E {x, y, z) = e^+^E (z) , H (x, y, z) = e i(k * x+k « y) H (z) , (38) 

which is a standard choice for the scattering problem of a plane electromagnetic wave incident 
on a plane-parallel stratified slab. Indeed, in such a case, due to the boundary conditions 
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()45j) . the tangential components (k x ,k y ) of the wave vector are the same for the incident, 
reflected and transmitted waves. The substitution (J3~%j) in Eq. (J3T|) allows to separate the 
tangential field components E x , E y , H x , H y into a closed system of four linear differential 
equations 
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d z ty (z) = i—M (z) ^ (z) , where \I> (z) 

c 



E x (z) 



Ey(z) 

H x (z) 
H y (z) 



(39) 



(40) 



The 4x4 matrix M (z) is referred to as the Maxwell operator. The reduced Maxwell 
equation for the four tangential field components \l/ (z) should be complemented with 
the following expressions for the normal components of the fields 

E z = {-n x H y + n y H x - e* xz E x - e* JZ E y ) e" 1 , 
H z = (n x E y - n y E x - [i xz EL x - V* yz H y ) 

where 

n x = ck x /u, n y = ck y /uj. (41) 

The expression for the Maxwell operator M (z) is very cumbersome, its explicit form can 
be found in The matrix elements of M (z) depend on the following parameters: the 
frequency u, the direction n of light incidence, the material tensors e (z) and fx (z). 

In a periodic layered medium 



M(z + L) =M(z) 



(42) 



where L is the stack period. For any given k x , k y and u, the system (|39p of four ordinary 
linear differential equations with periodic coefficients has four Bloch solutions 



* fcj (z + L) = e^ kt (z) 



1,2,3,4 



(43) 



where fcj, i = 1, 2, 3, 4 correspond to four solutions for k z for g 

k x , k y , u < — > {leu, k 2z , h z , hz} = {k$ z , k* 2z , kl z , k* Az } . (44) 

Real k z in (jUj) relate to propagating Bloch eigenmodes, while complex k z relate to the 
evanescent modes. In the case of propagating eigenmodes, the correspondence between k z 
and uj for fixed k x , k y is referred to as the axial dispersion relation, the concise form of which 
is given by Eqs. (JHJ) or ©. 
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The reduced Maxwell equations ([59]) in periodic layered media are analyzed and solved 
using the transfer matrix formalism, the detailed description of which can be found in the 
extensive literature on the S n bj ect (see, for example, SU HQ and re fe r enees therein). 



B. Boundary conditions 



The boundary conditions at the slab/vacuum interface reduce to the continuity require- 
ment for the tangential field components at z = 



(E, (x,y,0)) ± + (E R {x,y,0)) x 
[Hj (x,y,0)) ± + (H R (x,y,0)) ± 



(45) 



^Et (x,y,0)) ± , 
(H T (x,y,0)) ± , 

where the indices I, R and T denote the incident, reflected and transmitted waves, respec- 
tively Using representation (|5H)l . we can recast (j43j) in a compact form 



/ (0) + V R (0) = V T (0) 



(46) 



where 



R 



Ei,x 

E i,v 
Hi,x 

H i,y 

Er, x 

Er, v 
Hr, x 



(47) 



L H R,y 



Ei, x 

E I,y 

-E hx n x n y n; x - E ItV (1 
Ei,x{l-nl)n~ l + E lA 

Er, x 

E R, y 

E R)X n x n y n- 1 + E R>y (1 - n 2 x ) n' 1 
-E R)X (l - n 2 y ) n~ l - ER^nyn- 1 _ 

describe the incident and reflected waves, respectively. n x , n y , n z are the Cartesian compo- 
nents of the unit vector (fTU|) . 

Knowing the eigenmodes P5|) inside the slab and using the boundary conditions (f4T)|) 
we can express the amplitude and composition of the transmitted wave and reflected 
wave ^r, in terms of the amplitude and polarization of the incident wave This gives 
us the electromagnetic field distribution {z) inside the layered medium, as well as the 
transmittance and reflectance coefficients of the semi-infinite slab as functions of the incident 
wave polarization, the direction n of incidence, and the frequency u. 
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C. Energy flux, reflectance, transmittance 



The real-valued Poynting vector is denned by 

S (x, y, z) = ^3? [E* (x, y,z)xH {x, y, z)\ . (48) 
Plugging the representation (|3*5)l for E (x, y, z) and H (x, y, z) in Eq. (}4*%j) yields 

S(x,y,z) = S(z) = E* (z) x H (z) , (49) 

implying that none of the three Cartesian components of the energy density flux S depends 
on the tangential coordinates x and y. In addition, the energy conservation argument implies 
that the axial component S z of the energy flux does not depend on the coordinate z either 

S z (x, y, z) = S z = const, S x (x, y, z) = S x (z) , S y (x, y, z) = S y (z) . (50) 

This only apply to the case of a plane monochromatic wave incident on a lossless layered 
medium. The explicit expression for the z component of the energy flux ()49)1 is 

s z = \ [e*h v - e;h x + e x h; - e v h:} . (51) 

Let us turn to the scattering problem for semi-infinite slab. Let Si, Sr and St be the 
Poynting vectors of the incident, reflected and transmitted waves, respectively. The energy 
conservation imposes the following relation between the normal components of these three 
vectors 

'St) = [Si] + (S R ) . (52) 



Since the stack is presumably composed of lossless materials, the z component of the energy 
flux is independent of coordinates both inside and outside the stack. In particular, inside 
the slab we have 

at z > 0: {^Tj — const. (53) 
The transmittance (r) and the reflectance (p) of a lossless semi-infinite slab are defined as 



St) [Sr 
Si) ( Si 
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